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A mechanism of disappearance and formation of the Efimov levels of the 
helium 4 He trimer is studied when the force of the interatomic interaction is 
changed. It is shown that these levels arise from virtual levels which are in 
turn formed from (quasi) resonances settled on the real axis. The resonances 
including virtual levels are calculated by the method based on the solution of 
the boundary value problem, at complex energies, for the Faddeev differential 
equations describing the scattering processes (2 + 1 — > 2 + 1; 1 + 1 + 1). 
All the calculations are performed with the known interatomic Aziz He-He- 
potential HFD-B. A very strong repulsive component of this potential at short 
distances between helium atoms is approximated by a hard core. A special 
attention is paid to the substantiation of the method used for computing 
resonances and to the investigation of its applicability range. 



I. INTRODUCTION 

The 4 He three-atomic system is of considerable interest in various fields of physical 
chemistry and molecular physics. Studies of the helium dimer and trimer represent an 
important step towards understanding the properties of helium liquid drops, superfluidity in 
4 He films, and so on (see, for instance, Refs. [lr||). Besides, the helium trimer is probably a 
unique system where a direct manifestation of the Efimov effect [[| can be observed since the 
binding energy a of the 4 He dimer is extremely small (~ —1 mK 0-0) even in the molecular 
scale. For this reason, the helium trimer is certainly of interest for nuclear physicists, too. 
Moreover a theoretical study of the 4 He trimer is based just on the same methods of the 
theory of few-body systems that are used in solving three-body nuclear problems. 

From the standpoint of the general theory of few-body systems, the 4 He trimer belongs 
to three-body systems that are most difficult for a specific investigation, first, owing to its 
Efimov nature, and second, because it is necessary to take into account the practically hard 



core in the interatomic He - He - interaction B-ll|. At the same time the problem of three 
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helium atoms can be considered as an example of an ideal three-body quantum problem 
since the 4 He atoms are identical neutral bosons with zero spin and the analysis of this 
problem is complicated neither by separation of spin-isospin variables nor by taking into 
account the Coulomb interaction. 

There is a great number of experimental and theoretical studies of 4 He clusters. However, 
most of the theoretical investigations consist merely in computing the ground-state energies 

Besides, 



of clusters of that sort, mainly on the basis of variational methods [12-16 



the 



methods based on hyper spherical expansions of the Schrodinger and Faddeev equations 
fT7HT9| in the coordinate representation were used. Also, the Faddeev integral equations in 
the momentum representation were employed in Refs. |^O|,0 while the results of Ref. [f2^] are 
based on a direct solving the two-dimensional Faddeev differential equations in configuration 
space. From the experimental studies we would like to mention those of Refs. |5|-|7,23] where 
clusters consisting of a small number of noble gase atoms were investigated. 

Though much effort was undertaken for studying molecular clusters various problems re- 
lated to the 4 He trimer remained beyond the scope of thorough consideration. In particular, 
the elastic scattering phases of a helium atom on a helium dimer and breakup amplitudes (at 
ultralow energies) have been calculated only recently |2~4^[26f . These computations were pre- 
ceded only by the computation of characteristics of the He-He2 scattering at zero energy [2C 
and estimation of the recombination rate (1 + 1 + 1 — ^ 2 + 1) |27|| . 

As a matter of fact, we have already pointed out basic reasons for computations of excited 
states and scattering being especially difficult in the 4 He3 system. First, this is a low energy 
of the dimer which necessitates to consider very large domains in the configuration space 
with a characteristic size of hundreds of A. Second, a very strong repulsive component in the 
He-He interaction produces large errors in the standard approximation of the three-atomic 
Hamiltonian at short distances between atoms. The capacities of modern computers do not 
yet allow one to reach dimensions of grids that would remove both the above-mentioned 
reasons and would provide stable results with the use of the conventional methods. 

The present paper is a sequel of studies of the 4 He3 system undertaken in the pa- 
pers H24]-f26f within an approach that is capable, as we think, to resolve both the above- 
mentioned numerical problems. In these papers the repulsive component of the He-He 
interaction at short distances between atons is approximated by a hard core. This allows 
one to investigate the 4 He3 system within a mathematically rigorous method of solving a 
three-body problem in the Boundary-Condition Model developed in |28|,|29|]. An important 
advantage of such an approach that essentially diminishes computational errors is the ne- 
cessity to approximate, inside the core domains, only the Laplacian operator instead of the 
sum of this operator and a huge repulsive components of the He -He -potentials (see ||26|j). 
In |p4[-p6|l , such an approach has been successfully applied for calculating not only scattering 
but also binding energies of the ground and excited states of the helium trimer. Investigation 
made in fl24|-f2~6[| has shown that the method proposed in p8| , |29|| is well suited for performing 
three-body molecular computations in the case where repulsive components of interatomic 
interactions are of a hard core nature. 

There is a series of works |TB| , pT| , p6] showing that the excited state of the 4 He trimer is 
initiated indeed by the Efimov effect [Q. In these works the various versions of the Aziz 
4 He- 4 He potential were employed (HFDHE2 f§, HFD-B @, and LM2M2 0). However, 
the basic result of Refs. |T8|j21|j26[| on the excited state of the helium trimer is the same: this 
state disappears when the interatomic potential is multiplied by the "amplification factor" 
A of order 1.2. More precisely, if this potential is multiplied by the increasing factor A > 1 
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then the following effect is observed. First, the difference e^(A) — E\ (A) between the dimer 
energy e^(A) and the energy of the trimer excited state .E^(A) increases. Then the behavior 
of this difference radically changes and with further increase of A it monotonously decreases. 
At A ~ f .2 the level E^ disappears. It is just such a nonstandard behavior of the energy 
E^ 1 \\) as the coupling between helium atoms becomes more and more strengthening, points 
to the Efimov nature of the trimer excited state. And vice versa, when A slightly decreases 



(2) 

[no more than 2 %), the second excited state El appears in the trimer ||18l|21 



This paper is aimed at elucidating the fate of the trimer excited state upon its disap- 
pearance in the physical sheet when A > 1 and at studying the mechanism of arising of 
new excited states when A < I. As the interatomic He -He potential, we use the potential 
HFD-B We have established that for such He - He - interactions the trimer excited level 
E^p merges with the threshold at A ~ 1.18 and with further decreasing A it transforms 
into a virtual level of the first order (a simple real pole of the analytic continuation of the 
scattering matrix) lying in the unphysical energy sheet adjoining the physical sheet along 
the spectral interval between and the three-body threshold. We trace the position of this 
level for A increasing up to 1.5. Besides, we have found that the excited (Efimov) levels for 
A < 1 also originate from virtual levels of the first order that are formed in pairs. Before a 
pair of virtual levels appears, there occurs a fusion of a pair of conjugate resonances of the 
first order (simple complex poles of the analytic continuation of the scattering matrix in the 
unphysical sheet) resulting in the virtual level of the second order. 

As it will be clear from the further exposition (see Sect. |T|), the above-mentioned res- 
onances are not, generally speaking, genuine resonances of the 4 He3 trimer since they are 
situated outside of the energy domain for which we can rigorously prove the applicability of 
the method we are using for computing the resonances. We will call the resonances found 
outside the range of guaranteed applicability of the method the (quasi) resonances. 

The paper is organized as follows. 

In Sect. U, we describe the method of search for resonances in a three-body system on the 
basis of the Faddeev differential equations. The idea of the method consists in calculating 
the analytic continuation of the component Sq(z) (see formula (|6])) of the scattering matrix 
corresponding to the (2 + 1 — > 2 + 1) process, in the physical sheet with the use of these 
equations. A particular attention in this section is paid to the description of the parabolic 
domain on the physical sheet where one can analytically continue the function Sq(z) by 
numerical solving the coordinate space Faddeev partial equations. For the potentials we 
use, the three-body resonances (including virtual levels) lying in the unphysical sheet of 
energy z plane adjoining the physical sheet along the interval (e^, 0) are the roots of the 
function So(-2) in the physical sheet. We have earlier employed this method for computing 
resonances as roots of S (z) in the three-nucleon problem [pPf]. 

In Sect. fTJ, we first briefly describe the numerical method we use to solve the (2 + 1 — > 
2 + 1; 1 + 1 + 1) scattering problem for the 4 He3 system with going out into the domain of 
complex energies. Then we describe the results of our calculations. 

Some notation used throughout the paper is as follows: by C we denote the complex 
plane; y/z stands for the main branch of the function z 1 ^ 2 , Im^/z > for any z G C; the 
symbol IR+ is used for the quadrant x > 0, y > 0; by L 2 (M+) we understand the Hilbert space 
of complex-valued functions which are integrable on M. 2 ^ with the absolute value squared; 
the symbol ~z stands for the complex number conjugated to z. 
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II. METHOD FOR SEARCH OF RESONANCES IN A THREE BODY SYSTEM 
ON THE BASIS OF THE FADDEEV DIFFERENTIAL EQUATIONS 



A. Faddeev partial differential equations in the case of smooth potentials 

In this paper, we will consider the 4 He3 system in the state with the total angular 
momentum L = 0. 

First we consider the case where the interatomic interactions are described by conven- 
tional smooth potentials that include no hard-core component. In this case, the angular 
partial analysis reduces the initial Faddeev equation for three identical bosons to a system 
of coupled two-dimensional integro-differential equations |3l[ 

[H ,i - z\ Ffay) = -V(x)%(x,y) . (1) 

Here, x, y stand for the standard Jacobi variables, x > and y > 0, and 

d 2 d 2 ( 1 1 \ 

for the partial component of the kinetic energy operator. Functions from the domain of i?o,z 
are assumed to obey the boundary conditions 

Fi(x,y) \ x=0 = 0, F l (x,y)\ y=0 = 0. (3) 

which are quite standard when the expansions over bispherical basis are used. The potential 
V(x) is assumed to be central. In our paper, the energy z can get both real and complex 
values. At L = the partial angular momentum I corresponds both to the dimer and an 
additional atom. The momentum / can assume only even values, / = 0, 2, 4, . . . . 

The partial wave functions ^i(x,y) are expressed through the Faddeev partial compo- 
nents Fi(x,y) by the relations 

V l (x,y) = F l (x,y) + J2 d V h w {x,y, V ) F v {x' ,y') (4) 

v 

where 



, 1 2 3 2 v^3 . /3 2 1 2 v^3 

x = \h x -y^. y = y ^ ^ + ^ v +^ x yv- 



and — 1 < r] < 1. The explicit form of the functions hu> can be found, e.g., in Refs. |31|j32 
(see also Here we only deal with a finite number of equations ([!]), assuming that 

I < ^max where Z max is a certain fixed even number, Z max > 0. The condition < / < Z max is 
equivalent to the supposition that the potential V(x) only acts in the two-body states with 
I — 0,2, ... , / max - The spectrum of the Schrodinger operator for a system of three identical 
bosons with such a potential is denoted by cr 3B . 

Its is well known (see, e. g., Ref. that if the potential V(x) is smooth and decreasing 
as x — > oo together with its derivatives not slower than x~ 3 ~ £ , e > 0, then the asymptotic 
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conditions as p — > oo and/or y — > oo for the partial Faddeev components of the (2 + 1 — > 
2 + 1; 1 + 1 + 1) scattering wave functions^ for z = E + iO, E > 0, read 

Fi(x,y;z) = 5 lo ip d (x) {sm(y/z - e d y) + exp(iy/z - e d y) [& (z) +o(l)]} 

+ ^ W: ,,) + „)1, (5) 
VP 

We assume that the 4 He2 dimer has an only bound state with an energy e d , t d < 0, and 
wave function tp d (x). This function is assumed to be normalized so that for all x > values 

y 

of ijj d (x) are real. The notations p, p = \/ x 2 + y 2 , and 9, 9 = arctg — , are used for the 

hyperradius and hyperangle. The coefficient ao(z), 2 = -E + iO, for E > e d is the elastic 
scattering amplitude. The functions Ai(E + iO, 6*) provide us, at E > 0, the corresponding 
partial Faddeev breakup amplitudes. Note that for z = E + iO the correction terms o(l) in 
coefficients of outgoing waves exp(i v / ^ — cay), E > e d , and exp(iy / ip)/y / p, > 0, in (H) 
are of the form, respectively, o(y^ 1 ^ 2 ) and o(p -1 / 2 ). This property ensures uniqueness of the 
solution of the boundary value problem (|l]-||) for real scattering energies E > e d 

The (2 + 1— >-2 + 1) component of the s-wave partial scattering matrix for a system of 
three helium atoms is given for real z = E + iO, E > e d , by the expression 

S (z) = 1 + 2ia (^) (6) 

while the (2 + 1 — > 2 + 1) scattering phases read 

6o(p) = ^Im lnS (e d +p 2 + iO), p > 0, 
where p stands for the momentum conjugated to the Jacobi variable y. 



B. Holomorphy domains of the Faddeev components Fi{z) and scattering matrix So(z) 

Our goal is to study the analytic continuation of the scattering matrix So (z) into the 
complex plane (the physical sheet). As it follows from the results of Refs. |p3|j34}| , roots of 
the function S (z) in the physical sheet of energy z plane correspond to the location of the 
three-body resonances situated in the unphysical sheet connected with the physical sheet 
by crossing the spectral interval (e^, 0). This statement is a particular case of more general 
statements regarding the three-body resonances obtained in |3l|[34]] for the case of two-body 
potentials decreasing in the coordinate space not slower than exponentially. We assume that 
V(x) is just a potential which falls off exponentially and, thus, for all x > 

\V(x)\ <Cexp(-^), (7) 

with some positive C and p. For the sake of simplicity we even assume sometimes that V(x) 
is finite, i.e., V(x) = for x > r , r > 0. Looking ahead, we note that, in fact, in our 



1 Here we speak about the wave functions usually denoted by sign "(+)". The asymptotics of 
these functions in the total three-body configuration space R 6 contains, apart from the incident 
wave, only the so-called outgoing spherical waves (see, e.g., fl3l| ). 
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numerical computations of the 4 He3 system at complex energies we make a "cutoff' of the 
interatomic He - He - potential at a sufficiently large radius Tq. 

It is well known that different representations of the same holomorphic function (for 
instance, either by a series or by an integral) allow one to describe this function only in 
some parts of its Riemann surface. The description |$3|,|34| of the holomorphy domains for 



different truncations of the total three-body scattering matrix in the physical sheet was based 
on the use of the Faddeev integral equations in the momentum representation. In this paper, 
we make use of the Faddeev equations in the configuration space. Therefore it is necessary 
to perform an investigation, independent of |33| , |34|| , of domains in the physical sheet where 



we can analytically continue the Faddeev components Fi(x, y;z) and the amplitudes &o(z) 
and Ai(z, 6) just with the use of the configuration space techniques. 

Let us list briefly the main results of this investigation obtained by us for the 4 He3 
system under the assumption (|7|). To formulate these results we distinguish the following 
three domains in the complex plane C. 

1°. The domain n'*) where the Faddeev components Fi(x,y;z) (and, hence, the wave 
functions ^i(x, y; z)) can be analytically continued in z so that the differences 



$i(x,y;z) = Fi(x,y;z) - 5 w ip d (x) sm(y/z - e d y) (8) 

at z G n(*' \ o~3B turn out to be elements of L 2 (IR+). The domain n(*) is described by the 
inequality 



lmy/z -e d < min |— /i, V3VM j • (9) 

For fixed x, y the functions <&i(x, y;z) are continuous in z up to the rims of the cut along 
the continuous spectrum [e^, +oo). 

2°. The domain 11^ where both the elastic scattering amplitude ao(z) and the Faddeev 
breakup amplitudes Ai(z, 6) can be analytically continued in z, z £ a 3B , provided that the 
functions Fi(x,y;z) obey the asymptotic formulas (|5]). This domain is described by the 
inequalities 

I y^3 

Im yfz + - Im yjz - e d < — VM > ( 10 ) 



2 

Im \fz + Im \J z — e d < -—-/i . (If) 




3°. And finally, we distinguish the domain U^ s \ most interesting for us, where the 
analytic continuation in z, z ^ a^B, can be only done for the elastic scattering amplitude 
& (z) (and consequently, for the scattering matrix S (z)); the analytic continuabilty of the 
amplitudes Ai(z,9) in the whole domain is not required. The set is a geometric 
locus of points obeying the inequality 




Im y/z - e d < min < —= V |e d |, —pr-fi > ■ (12) 
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For the domains 11^, 11^, and U^ s \ the following chain of inclusions 

U (A) c n (S) c nW _ 

is valid. 

Note that the type @ or (|T2|) condition, 

Im v^i - a < A a G 1, 6 > 0, (13) 
is equivalent to the inequality 

Rez > a - b + — (Imz) 2 . (14) 

Therefore, for jj < 2y/\e,i\ the set is the domain bounded by the parabola 

3 1 

Rez > -\e d \ - -a 2 + — (Imz) 2 . (15) 

4 i\i A 

For \i > 2y/\ed\ this set coincides with the domain 

Rez> -4|e d | + — l —(\mz) 2 . (16) 

Analogously, if fi < then the domain is described by the inequality (|15|); whereas 

for /i > | yfidj by the inequality 

Re^>-^|e d | + -^(Imz) 2 . (17) 

As to the curves bounding the domains (|l0|) and (|TTD , we only notice that their order with 
respect to the variables Re z and Imz is higher than the second order. It is easy to check 
that each of these curves is connected, symmetric with respect to the Re z axis and crosses 
the latter only once. For the first curve this intersection occurs at z = |e^, the slope angle 



= As 

lmz=±0 



of the tangent at the point of intersection being independent of e d , i/h 

Rez — > +oo, the boundaries ( |TTj ) and flTT|) are asymptotically approximated by the type ( |T4"D 
parabolas with coefficients a and 6 which can be computed explicitly. 

To prove the assertion 1° concerning the domain II^) we note that the functions 
&i(x, y; z) given by the formulas (H) satisfy the equations 

r+i 

[H 0t i+V(x) - z]$i{x,y;z)+V{x)y] / d-qh w {x,y,-q) $i>{x',y',z) = xi{x,y;z) (18) 

v J - 1 

where 

-+1 



Xi{x, y\ z) = -V(x) j drj h m (x, y, r])ip d (x') sin(y/z - e d y') 
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Obviously, for z G H^' the functions (a:, j/; z) fall off exponentially as p — ► oo. Moreover, 
for all the directions < 9 < tt/2 the uniform estimate 



Xi(x,y;z) = 0(exp(-ap) 



is valid with a = min — Im y/z — e^, — I ^ m V 2 — e czj • Consequently, if the 

condition (|9]) holds then the inhomogeneous terms Xi{ x )V'-> z ) considered as functions of the 
variables x and y at fixed z, are elements of L 2 (IR+). At the same time, the vectors Xl( z ) 
turn out to be holomorphic functions of z G IT*) with respect to the L 2 (R+) norm. 

In the problem under consideration, the spectrum of the Faddeev matrix operator defined 
by the l.h.s. of Eqs. ( |T%D and by the boundary conditions (Q) in the Hilbert space constituted 
of the vectors <3> = ($ , $2, • • • , $i miX ), $i G L 2 (M.+), coincides with the spectrum a 3B of the 
corresponding three-boson Schrodinger operator with two-body potentials V(x) only acting 
in the states with I — 0,2, ... , Z max . This means that for any energy z G lying outside 
of the spectrum a 3B , the inhomogeneous system ([TJ) is uniquely solvable in the class of 
the functions ®i(z) G L 2 (M+), I — 0,2, ... , Z max . Since outside of the set a 3B the resolvent 
of the Faddeev operator is a holomorphic operator-valued function of the variable z, each 
of the components $i(z) of the solution of Eqs. QT8D also is a holomorphic function of 
z G \ a 3B . The bound-state energies of the three-boson system under consideration 
turn out to be poles of the first order for $i(z). Thus, the Faddeev partial components 
Fi(x,y; z) admit the analytic continuation in z in the form into the domain IT*) \ <j 3 b- 

The proof of the assertions of 2° and 3° regarding the domains and IT-^ is rather 
cumbersome. This is why we here only outline its main steps. Note that the proof is based 
on the integral equations method and it is quite standard (see, for instance, Ref. ||31|| , Chap- 
ter V). First, the equations fll8|) are rewritten in the form of the Faddeev partial integral 
equations. To do this, it suffices to reverse the operators if ,z + V — z in (|18D. Since the 
variables x and y in H 0y i are separated, the kernels {X, X'; z), X = {x, y}, X' = {x', y'}, 
of the respective resolvents Rv\z) = (H i + V — z)' 1 are explicitly expressed in terms of 
the two-body problem. Analytic properties in the variable z and coordinate asymptotics of 
the kernels R^(X, X';z) are well known (see Ref. ||31|| , Chapters IV and V). Iterations first 



"improve" and then stabilize the asymptotic properties of the iterated kernels and inhomo- 
geneous terms of the Faddeev equations. (In the case under consideration, this stabilization 
requires only three iterations.) Further, it turns out that, for z G" [e^, +00), the iterated ker- 
nels are represented by sums of exponentially decreasing terms admitting, in certain domains 
of the configuration space, an explicit asymptotic factorization with respect to X and X'. 
Since we are working in the domain where Xi{ z ) £ ^2^+); the corresponding asymptotic 
factors of these terms, along with the asymptotics of the iterations of the inhomogeneous 
term, determine the coordinate asymptotics of the functions $>i(x,y;z). Therefore, finally 
we are able to determine the geometric locus of the points z in the complex plane for which 
there exists a (non-empty) set in the configuration space such that the leading term of the 
coordinate asymptotics of the function &i(x,y; z) in this set represents a term of the form 
&o(z)%l)d(x) exp(i v / 2 — €dy), and thus, for these z the scattering matrix S (z) is well defined. 
This geometrical locus is just the domain iT 5 ). In this domain, as y —* 00 and/or p — > 00, 
the functions $>i(x, y,z) admit the asymptotic representation 

$i(x,y;z) = 5i ipd{x) {exp(iy/z - e d y) [a (z) + o(l)] + f (y; z)} (19) 
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+ eMi f p) [Mz, 6) + (1)] + F hl (x, y; z) 

VP 

with 

fofaz) = O(e- ao ^ y ) and F u (x,y;z) = (e~ ai ^ p ) (20) 

y— >oo ' p— >oo 

where 

«o(^) = — VN - Im^-fd, (21) 

ai(z) = min |«o(-2) , At - Im - ed , Im-\/i| • (22) 

In a parabolic neighborhood of the y-axis, the functions F lt i(x, y; z) are also subjected to 
the asymptotic estimates 

F hl (x,y;z) = O (exp(-a (z)y)j (23) 
y — > oo v 7 

x < 

where z/ is an arbitrary fixed number smaller than unity, v < 1. 

As to the domain n^, the leading asymptotic term of each of the functions Fij(x, y; z) 
for z G IT > is a spherical wave exp(iy / ip)/y / p with the amplitude A\j(Q) being a dif- 
ferentiable function of the angle 9. Therefore, for z G IT^ the term Fi i /(x, y; z) in the 
r.h.s. of the formula flT9|) can be added to the asymptotic term with a spherical wave 
preceding F^. In the domain Il( 5 ), and hence, in a narrower domain Il( A ) the condition 

f (y; z) = o[ exp(i^z — e d y) ) holds. Consequently, for z G 11^, the Faddeev components 

y^oo \ J 

Fi(x,y; z) do obey the standard asymptotic conditions like (||). 

Therefore, for any v < 1 the dominant term of the asymptotics of the function <&o(x, V'-, z )i 
z G IiS s ' \ <73£, in the domain x < y u reads as a (z)ipd(x) exp(iy/z — e^y) as y — > 00. 
This means that, for z G IT^ \ a^B, it is always possible by solving the equations ([[]) to 
separate explicitly the elastic scattering amplitude &q(z) and, thus, to construct the analytic 
continuation of the scattering matrix S (^). 

Outside of the domain iF 5 ) the numerical construction of So (z) by solving the Faddeev 
differential equations is, in general, impossible since for x < y v and v < 1 both functions 
io{y\z) and F 1>0 (x, y; z), z G" IT^, include terms decreasing slower than exp(iy/z — y) as 
y -> 00. 

C. The partial Faddeev differential equations in the case of potentials with hard core 

In the case of potentials with hard core, the partial Faddeev differential equations for a 
system of three identical bosons at L = acquire the form 

W»- MM = {- Vlxy H?* v) ' HI (24) 
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where c, c > 0, is the core size. The partial wave functions *&i(x, y) are expressed via 
Faddeev partial components Fi(x,y) by the formulas (|j). The components Fi(x,y) satisfy 
the standard boundary conditions (S). The two-body central potential V(x) acts only beyond 
the core domain, i. e. only where x > c. We assume as before that V(x) falls off not slower 
than exponentially as x —>■ oo and, hence, it satisfies the condition (0) for some C > and 
H > 0. 

A main difference between the model with hard core and those with smooth potentials 
is that the functions Fi(x,y) in this model satisfy the auxiliary boundary conditions 

F t (c,y) + J2 f d V h w (c,y, V )F l ,(x',y') = (25) 
i' 

requiring that the wave functions ^i(x, y) vanish on the boundary x = c of the core domain. 
It can be shown that in fact the conditions (|25"D force the wave functions (|J) to vanish also 
inside the core domain at all energies except for a certain countable set of real values of z 



(see Ref. [BE] and references therein). 

Asymptotic conditions for the partial Faddeev components Fi(x, y;z) of the (2 + 1 — > 
2 + 1; 1 + 1 + 1) scattering wave functions as p — > oo and/or y — > oo are again of the form 
(H). The only difference is that the dimer wave function ip d {x) is considered as zero in the 
core domain, i. e. for < x < c. 



In the hard-core model, all the assertions of Sect. II B regarding the holomorphy domains 



of the functions &i(x, y; z) and the scattering matrix So (z) in z still hold true. 



D. Resonances and virtual levels as roots of the scattering matrix So(z) in the 

physical sheet 

We have already noticed that the roots of So (z) in the physical sheet of energy z plane 
correspond to the location of the three-body resonances in the unphysical sheet adjoining 
the physical sheet along the spectral interval (e<j,0). In the case under consideration, this 
statement is an immediate consequence of the unitarity of the scattering matrix Sq(z) for 
z = E + iO, e d < E < 0, 

S (£ + i0)S (£ + i0) = 1. (26) 

Indeed, as we have established, the functions $>i(x, y;z) are holomorphic functions of z G 
\ o"3B- Since the boundary value problem (jl]-|5D is uniquely solvable, one easily verifies 
that the boundary values &i(x, y \ E + iO) and $;(x, y;E — iO) for these functions on the rims 
of the cut along [e d , +oo) are related to each other as 

^(x,y;E + i0) = y- E - iO) (27) 

since, on the one hand, their asymptotics (O) as y — ► oo and/or p — ► oo has the same 
structure and, on the other hand, 

Xi(x, y,E + iO) = -xi(x, y,E- iO) = ~xi(x, y;E - iO) , 

since 

sm(\/E - e d + i0y) = -sin( v / E - e d - iOy) = -sin(y / E - e d - iOy) . 
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Consequently, 

& {E + iO) = -& {E - iO) (28) 

and 

S (£ + iO) = S (£-iO), E>e d . (29) 

Therefore, it follows from Eq. ( Pop that for < E < 

S (£ + iO) = [So^-iO)]- 1 and S (E - iO) = [S (E + iO)]" 1 . 

This means that the function So(-E + iO) is continued through the cut [e<j, 0] into the domain 
lm.z < ). In a similar manner, Sq(E — iO) is continued into the domain Imz > 0, 

again as Sq 1 (z). All this signifies that the scattering matrix Sq(z) admits analytic continua- 
tion at least into the domain II^ of the unphysical energy sheet connected with the physical 
sheet by crossing the interval [e^, 0], the value of the continued function Sq(z) at z G IT^ in 
the unphysical sheet coinciding with the value of Sq (z) at the same z but in the physical 
sheet. 

Recall that those points z on unphysical sheets are called resonances where the analyt- 
ically continued scattering matrix possesses poles. The resonances z with zero imaginary 
part Im z = and Re z < €d are called the virtual levels. 

Thus, we have here presented a simple proof of the fact that the resonances including the 
virtual levels corresponding to poles of the analytic continuation of the scattering matrix 
Sq(z) in the unphysical sheet connected with the physical one by crossing the spectral interval 
[e<2, 0] are the roots of this matrix in the physical sheet. At the same time, the poles of the 
function Sg(z) in the physical sheet correspond to bound states of the three-boson system 
under consideration. 

Concluding the subsection, we note that it follows from Eq. (|28|) that &q(z) = — a (z) 
and, hence, 

S^) = So(z) (30) 

for any z G 11^. This means that the roots of the function Sq(z) are situated symmetrically 
with respect to the real axis. 



III. NUMERICAL METHOD AND RESULTS OF COMPUTATIONS 



In the present work we make use of the Faddeev equations (24) considered together with 
the boundary conditions (§), @ and (p5|) to calculate the values of the 4 He3 scattering 
matrix Sq(z) in the physical sheet. We search for the resonances including the virtual levels 
as roots of So (z) and for the bound-state energies as positions of poles of So (z). All the 
results presented below are obtained for the case / max = 0. 

In all our calculations, H 2 /m = 12.12 KA 2 . As the interatomic He -He- interaction we 
employed the widely used semiempirical potential HFD-B constructed by R. A. Aziz and 
co-workers 0. This potential is of the form 



Vhfd-b(z) = e { ^exp(-aC + P( 2 ) - 



Ce Cg Cio 



no 



(31) 
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where ( = x/r m . The function F(Q reads 

F(c) = { exp [- (D/C - l)] 2 , ifC<^ 
^ \l, i£C>D. 

For completeness the parameters of the potential HFD-B are given in Table I. 

The value of the parameter c (the core "diameter" of particles) is chosen to be so small 
that its further decrease does not appreciably influence the dimer binding energy and the 



energy of the trimer ground state Ef\ Unlike papers p4j-j26], where c was taken to be equal 



0.7 A, now we take c = 1.3 A. We have found that such a value of c provides at least six 
reliable figures of and three figures 

Since the statements of Sect. [II] are valid, generally speaking, only for the potentials 
decreasing not slower than exponentially, we cut off the potential HFD-B setting V(x) = 
for x > r . We have established that this cutoff for tq > 95 A provides the same values of 
e<2 (ed = — 1.68541 mK), (E^ = — 0.096 K) and phases 8q(p) which were obtained in 
our earlier calculations performed with the potential HFD-B. Comparison of these 

results with results of other researchers can be found in Refs. [p4| -|26]] . In all the calculations 
of the present work we take r = 100 A. Note that if the formulas from Sect. [H] including 
the parameter /i are used for finite potentials, one should set fi = +oo. 

Before making numerical approximation of the system of equations @, (p4|), (|25| ) at 
l max = we rewrite it in terms of a new unknown function <&o(x, y; z) that is expressed via 
the Faddeev component F (x,y;z) by the relation (g). Note that for z G U.^> \ ct 3 b the 
function $o(x, y;z) is square integrable in x, y (see Sect. |11 B| ) . Therefore, this function is 
uniquely determined by the asymptotic condition 

%(x,y;z) ^0 (32) 

p— >oo 

that can be easily approximated and programmed. One could, for instance, require 

§o(x,y; z) = at a sufficiently large p max and look for a numerical solution 

y :r 2 +j/ 2 =p max 



of the system (^), (p4]), (p5|) satisfying this condition. Further, for z G H^ s >, one could, 
going sufficiently far from p max into the domain of smaller (but nevertheless, providing the 
asymptotics (|19D) values of p, separate the elastic scattering amplitude ao(z), putting, e.g., 
ao(z) ~ &o{x, y; z) exp(— i^/z — ejy), where the value of x corresponds to the maximum 
of the function ipd{x). Such an approach is, however, not effective in view of a relatively 
slow decrease of the exponentials exp(— yjedj x) and exp(— Im y/z — t^y) as well as of the 
function exp(— Im \fz p) in the energy domain of interest for us in n( s ). For a proper ap- 
proximation of the condition (|32]), very large values of p max are to be taken. This is just 
a reason why one should take into account the asymptotics of the function $ (x,y;z) as 
x — ► oo and/or y — > oo. Though the asymptotic formula fl5|) only holds for z G we 
employ it also for z G \ U^ A \ Indeed, when z G \ U^ A \ the leading term of the 
asymptotics of §o(x,y; z) as y — » oo and x < y u , v < 1, is given by the same expression 
a (z) exp(i^/z — ejy) (see Sect. |11B| ) as in Eq. (|5]). Outside of the parabola x < y u , it suf- 
fices to require the condition ([$2]) to be satisfied. The presence, in Eq. (||), of the spherical 
wave exp(iy / i p)/y / p does not contradict this requirement. Therefore, the use of asymptotic 
condition (|j|) is justified even if z G n*- 5 -* \ H^ A \ 

A detailed description of the numerical method we use is presented in Ref. P5[ . Here 
we only mention main steps of the computational scheme |25| helpful for understanding our 
results. 
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When solving the boundary- value problem ([3D, @, ([24]), (p5|) written in terms of the 
function <3>o, we carry out its finite-difference approximation in polar coordinates p and 9. 
The grid is chosen in such a way that the points of intersection of arcs p — i — 1,2, ... , N p , 
and rays 9 — 9j, j — 1, 2, . . . , Ng, with the line x = c turn out automatically to be its knots. 
The pi points are chosen according to the formulas 



i = l,2,...,N^, 



c^ + yl i = l,2,...,N„-N^, 



where iVi^ stands for the number of arcs inside the core domain and 



Vi = f{n)sJp 2 Np - c 2 , r, 



' ~ N n - « ' 



The nonlinear monotonously increasing function /(r), < r < 1, satisfying the conditions 
/(0) = and /(l) = 1 is chosen in the form 



/(t) 



a i" , t G [0, r ] 
air + r v , r G (r , 1] 



The values of a®, Uq > 0, and cti, cti > 0, are determined via r and z/ from the continuity 
condition for /(r) and its derivative at the point To- As a rule, we took values of r within 0.1 
and 0.2. The value of the power v depends on the cutoff radius p max = pn p = 50 — 4100 A 
its range being within 2 and 4 in our calculations. 

The knots 9j at j — 1, 2, . . . , N p — Nc P ^ are taken according to 9j = arctg(?/j / c) . The rest 
knots 9j, j = N p — Nc P ^ + 1, . . . , Ng, are chosen equidistantly. Such a choice of the grid is 
prescribed by the need to have a higher density of points in the domain where the functions 
y; z) are most rapidly changing, i.e. for small values of p and/or x and lower in the 
asymptotic domain. In this work, we used the grids of dimension Nq = N p =600 — 1000. 
The number of the last arc knots in 9 lying in the core domain was usually equal to Nc P ^ = 5. 



The finite-difference approximation of the integro-differential equations ( p4"| ) and bound- 
ary conditions (|3|) , (^5|) for Z max = reduces the problem to a system of N p N g linear algebraic 
equations. The finite-difference equations corresponding to the arc i = N p include initially 
the values of the unknown function <3>(x, y; z) from the arc % — N p + 1. To eliminate them, 
we express these values through the values of y; z) on the arcs i = N p and i = N p — 1 
by using the asymptotic formula (H), just in the manner described in the concluding part 



of Appendix A of Ref. [26]. In [26], this approach was only used for computing the energies 



of bound states. Now we extend it also on the scattering problem. (Note that the formulas 



(A10) and (All) in |26| related to the described approach contain misprints. The values 
Cjj p in these formulas should be replaced with inverse values 1/Cjj .) The matrix of the 
resultant system of equations has a block-three-diagonal form (see Ref. f26|, Appendix A). 
Every block has the dimension Ng x Ng and consists of the coefficients standing at unknown 
values of the function y;z) in the grid knots belonging to a certain arc p = p^. The 
main diagonal of the matrix consists of N p such blocks. 

In contrast to P^j-pS] , in the present paper we solve the block-three-diagonal algebraic 
system on the basis of the matrix sweep method. This allows us to dispense with writing the 
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system matrix on the hard drive and to carry out all the operations related to its inversion 
immediately in RAM. Besides, the matrix sweep method reduces almost by one order the 
computer time required for computations on the grids of the same dimensions as in p3-GB| . 

We searched for the resonances (roots of the function Sq(z) on the physical sheet) and 
bound-state energies (roots of the function Sq l (z) for real z < €4) of the helium trimer by 
using the complex version of the secant method. Within this method, the approximation z n 
to a root of a holomorphic function f(z) is constructed from the two previous approximations 
z n ^i and z n _ 2 according to the formula 

_ _ f(z n -l){Zn-l ~ Zn-2) 

U ~ l f(Zn-l) ~ f(Zn-2)) 

As the relationship (^) implies the symmetry of properties of the scattering matrix Sq(z) 
with respect to the real axis, we performed all the calculations for So (2) only for Imz > 
(except the tests of the code). We start with a study of graph surfaces of the real and imag- 
inary parts of the scattering matrix So (z) in the domain of its holomorphy \ a^B ■ The 
root lines of the functions Re S (z) and Im Sq(^) obtained in the case of the grid parameters 
Ng = Np = 600 and p max = 600 A are depicted in Fig. 1. Both resonances (roots of Sq(z)) 
and bound-state energies (poles of So (z)) of the 4 He trimer are associated with the intersec- 
tion points of the curves Re So (z) = and Im So (z) = 0. When the roots or poles are simple, 
these curves intersect each other at the right angle. Note that for real z < e& the function 
So (2) is real and, thus, ImSo(^) = 0. In Fig. 1, along with the root lines we also plot the 
boundaries of the domains U^ s \ Il( A ) and IT*). One can observe that a "good" domain IT^ 
includes none of the points of intersection of the root lines ReSo(z) = and ImSo(z) = 0. 
Nevertheless, as we will see below, the going beyond this domain is of an interest, even 
though the asymptotic formula (0) is not valid for z G C\n^ and the function So (z) calcu- 
lated there cannot be interpreted as the scattering matrix. The caption for Fig. 1 points out 
positions of the four "resonances", the roots of So(-z), found immediately beyond the bound- 
ary of the domain n^. As one could expect, the values of the function So (2) at z E C\n( 5 ) 
and positions of its roots in C\Il( 5 ) turn out to be unstable and strongly depend on the value 
of the cutoff radius p max , whereas the dependence on the number of knots is weak. In par- 
ticular, for p max = 400 A, a (quasi) resonance, closest to the real axis, is situated at the point 
(-1.95 + i 1.81) mK, if Ng = N p = 300, at the point (-1.90 + i 1.85) mK, if Ng = N p = 520, 
and at the point (—1.89 + i 1.86) mK if Ng = N p = 800. The same (quasi) resonance in Fig. 1 
(calculated for p max = 600 A) is situated at the point (-2.34 + i 0.97) mK. If N e = N p = 600 
is fixed, the increase of p max up to 800 A shifts this point to the point (—2.44 + i 0.65) mK. 

All the aforesaid regarding the instability of the function So (z) values and positions of its 
roots beyond the domain IT^ bears no relation to its pole at the point z = E^ = —2.46 mK, 
corresponding to a trimer excited-state energy, even though this energy does not belong 
to IT 5 ). The point is that the position of the pole of So (2) is only determined by the 
position of the root of the determinant of the linear algebraic system we solve, whereas the 
inhomogeneous term of the system plays no role. Therefore, the search for the poles of the 
grid function So (z) is equivalent to the search for the binding energies of the trimer. The 
grids we have used turn out to be quite sufficient for this purpose. The convergence of 
our results for E^ with respect to the parameters Ng, N p , p max and their accuracy can be 
judged from the values of the difference — obtained with different grids and shown 
in Table II. 
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We would like to stress that we do not consider the roots of function Sq(z) drawn in 
Fig. 1 as genuine resonances for the He3 system since they are situated beyond the domain 
where the applicability of our method is proved. We should rather consider them as 
artifacts of the method. However it is remarkable that the "true" (i.e., getting inside 11^) 
virtual levels and then the energies of the excited (Efimov) states appear just due to these 
(quasi)resonances when the potential V(x) is weakened. This is the object of our further 
consideration. 

Following |18| , |21|j26[| , instead of the initial potential V(x) = Vhfd-bOe), we will consider 
the potentials 

V(x) = A • Vhfd-b(V)- 

To establish the mechanism of formation of new excited states in the 4 He trimer, we have 
first calculated the scattering matrix Sq(z) for A < 1. In Table III for some values of A from 
the interval between 0.995 and 0.975, we present the positions of roots and poles of S (z), 
we have obtained at real z < e^(A). We have found that for a value of A slightly smaller than 
0.9885, the (quasi)resonance closest to the real axis (see Fig. 1) gets on it and transforms into 
a virtual level (the root of So(z)) of the second order corresponding to the energy value where 
the graph of S (z), z G M., z < e<f, is tangent to the axis z. This virtual level is preceded by 
the (quasi)resonances z = (-1.04 + i 0.11) mK (z/\e d \ = -1.58 + i 0.168) for A = 0.989 and 
z = (-0.99 + i0.04)mK (z/\e d \ = -1.59 + i 0.064) for A = 0.9885. The originating virtual 
level is of the second order since simultaneously with the root of the function So (z), also 
the conjugate root of this function gets on the real axis. With a subsequent decrease of A 
the virtual level of the second order splits into a pair of the virtual levels E^* and E^ * , 
fi( 2 )* < ^( 2 )** Q £ ^ e first order which move in opposite directions. A characteristic behavior 
of the scattering matrix So (z) when resonances transform into virtual levels is shown in 
Fig. 2. The virtual level E^** moves towards the threshold e d and "collides" with it at 
A < 0.98. For A = 0.975 the function So(z) has no longer the root corresponding to E\ 
Instead of the root, it acquires a new pole corresponding to the second excited state of the 
trimer with the energy E^. Note that though the virtual levels E^* and E^** appear 
beyond the domain H^ s \ already at A = 0.985 the point E^* turns out to be inside this 
domain. Therefore, it should be considered as a "true" virtual level of the trimer. We expect 
that the subsequent Efimov levels originate from the virtual levels just according to the same 
scheme as the level E\ does. 

The other purpose of the present investigation is to determine the mechanism of disap- 
pearance of the excited state of the helium trimer when the two-body interactions become 
stronger owing to the increasing coupling constant A > 1. It turned out that this disappear- 
ance proceeds just according to the scheme of the formation of new excited states; only the 
order of occurring events is inverse. 

The results of our computations of the energy E^ when A changes from 1.05 to 1.17 are 
given in Table IV. In the interval between A = 1.17 and A = 1.18 there occurs a "jump" of 
the level Ef 1 on the unphysical sheet and it transforms from the pole of the function So(-z) 
into its root, E\ ' , corresponding to the trimer virtual level. The results of calculation of 
this virtual level where A changes from 1.18 to 1.5 are presented in Table V. For all the 
values of A presented in Tables IV and V, the dimer possesses an only bound state. We have 
found that the first excited state of the dimer appears only at A = 6.81. 

Note that in the case of finite potentials the geometric characteristics of the domain IT 5 ) 
where the function Sq(z) can be calculated reliably, are only determined by the value of 
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|ed(A)| (see formula (§) for p = +00). When |e<j(A)| increases, the domain IP 5 ) is enlarged. 
It is easy to check that the energies of the excited state level (A) and of the virtual level 
E^*(X) given in Tables IV and V belong to the corresponding domains n( s )(A). For A > 1, 



this results in a weak dependence of the calculated values of .E^(A) and E^'*(X) on the 
parameters N e , N p and (this is especially important) on the parameter p max . 

In essential, we chose the values of the cutoff hyperradius p max given in Tables III-V 
from the scaling considerations. As a matter of fact, we took the value of p max following the 
formula 



-.(I)*/ 



Pmax (A) 



Ci 



(33) 



where the "constant" C\ 



A=l 



corresponds to an appropriate choice of p max at 



400 



A=l 



A = 1. It has been established in [p4|-|26| that such a choice is ensured if p, 

600 A. In determining the values of p max (A), indicated in Table III, we followed the for 
mula (|33"D literally. As the "constant" C\, we took its value corresponding to the base 

value of p max =600 A. The values of p ma x(A) presented in Tables IV and V correspond 

in the interval within 400 and 800 A. All the results presented in 



A=l 

to the choice of p r 

Tables III-V have been obtained with the grids parameters Ng = N p = 600. 



A=l 
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I. The parameters for the HFD-B 4 He— 4 He potential. 



e(K) 

m (A) 

A 

a 

(3 
C 6 
Cs 
Cio 
D 



10.948 

2.963 

184431.01 

10.43329537 

-2.27965105 

1.36745214 

0.42123807 

0.17473318 

1.4826 



II. Dependence of the difference ed — E\ (mK) between the dimer energy and the trimer 
excited state energy E ( P on the grid parameters. The values of p max are in A. 



N 6 ,N p (Ng = Np) 

Pmax 


600 


800 


1000 


400 


0.7752 


0.7661 


0.7625 


600 


0.7809 


0.7695 


0.7649 


800 


0.7852 


0.7723 


0.7669 



III. The dimer binding energy e<z and the differences ed~E t , e^ — e\' ed — E t and ed — E\ 



(all in mK) between this energy and the trimer exited-state energies E\ , E t ' and the virtual-state 
energies E t , E\ ' depending on factor A. 



A 


£d 






e d ~ E?>* 


e d ~ E t 


Pmax (A) 


0.995 


-1.160 


0.710 








723 


0.990 


-0.732 


0.622 








910 


0.9875 


-0.555 


0.573 


0.473 


0.222 




1046 


0.985 


-0.402 


0.518 


0.4925 


0.097 




1228 


0.980 


-0.170 


0.39616 


0.39562 


0.009435 




1890 


0.975 


-0.036 


0.2593674545 


0.2593674502 




0.00156 


4099 



IV. Dependence of the dimer energy and the difference — E\ between this energy and 



the trimer exited-state energy E). on the factor A. 



A 


e d (mK) 


e d - E { 1 ] (mK) 


Pmax (A) 


1.05 


-12.244 


0.873 


300 


1.10 


-32.222 


0.450 


200 


1.15 


-61.280 


0.078 


150 


1.16 


-68.150 


0.028 


120 


1.17 


-75.367 


0.006 


120 
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( 1 )* 

V. Dependence of the dimer energy and the difference — E t between this energy and 



the trimer virtual-state energy E\ on the factor A. 



A 


e ( i (mK) 


e d - Ej 1] * (mK) 


//ni 'A y ( -Ti- ) 

1 llldA \ / 


1.18 


-82.927 


0.001 


110 


1.19 


-90.829 


0.016 


110 


1.20 


-99.068 


0.057 


100 


1.25 


-145.240 


0.588 


85 


1.30 


-199.457 


1.831 


70 


1.35 


-261.393 


3.602 


70 


1.40 


-330.737 


6.104 


55 


1.50 


-490.479 


12.276 


50 
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He z (mK) 

. 1. Root locus curves of the real and imaginary parts of the scattering matrix So (z). The solid 
lines correspond to Re So(^) = 0, while the tiny dashed lines, to ImSo(z) = 0. The Numbers 1, 2, 3 
denote the boundaries of the domains IlW, Tl^ and respectively. Complex roots of the func- 
tion So(z) are represented by the crossing points of the curves ReSo(z) = and ImSo(-z) = and 
are located at (-2.34 + i 0.96) mK, (-0.59 + i 2.67) mK, (2.51 + i 4.34) mK and (6.92 + i 6.10) mK. 
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